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Abstract 

We calculate the asymptotic behaviour of correlation functions as a function of the microscopic 
parameters for a Bose-Fermi mixture with repulsive interaction in one dimension. For two cases, 
namely polarized and unpolarized fermions the singularities of the momentum distribution func- 
tions are characterized as a function of the coupling constant and the relative density of bosons. 
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In the past years, the advances in cooling and trapping of atomic gases have opened 

the possibility to realize quasi one-dimensional (ID) systems with tunable strength of the 

interactions in optical lattices. This gives rise to new opportunities for the investigation of 

the striking phenomena appearing in correlated systems as a consequence of the enhanced 

quantum fluctuations in reduced spatial dimensions. The observable signatures of these 

phenomena are encoded in the correlation functions of the system such as the momentum 

distribution function which can be measured directly in time-of flight experiments or using 
rn 

Bragg spectroscopy A setup for the measurement of various density correlation functions 
has been proposed for the identification of dominant correlations in the atomic gas 
Theoretical studies of correlation functions in cold atomic gases have been performed both 
using analytical methods, e.g. bosonization p combined with exact results from integrable 
models such as the Bose-gas with repulsive 5-function interaction 4], and numerically. New 
correlation effects appear when the particles considered have internal degrees of freedom, 
n cold gases containing different constituent atoms Bose- Fermi mixtures can be realized 
sj]. Extensive theoretical results exist for ID Fermi gases due to their equivalence with 
the Tomonaga-Luttinger (TL) liquids realized by correlated electrons in ID lattices 0,0]. 
Only recently, theoretical investigations have been extended to Bose-Fermi mixtures: some 
correlation functions have been calculated numerically in the strong coupling limit [?], 0] 
where the problem simplifies due to the factorization of the many-particle wave function 
(see e.g. j2|). For analytical results on these systems one has to go beyond mean- field 
approximations and use methods which can capture the strong quantum fluctuations in ID 
systems. The phase diagram and certain correlation functions of atomic mixtures have been 
studied in the Luttinger liquid picture [lol.lll|. Without further input, however, these results 
are limited to the weakly interacting regime since the TL parameters which determine the 
low-energy theory cannot easily be related to the microscopic parameters describing the 
underlying gas. Therefore, instabilities predicted within this approach may not appear in a 
specific realizations jl, 12, 13^ . 

In this paper we establish the relation between the TL and the microscopic parameters 
for an integrable Bose-Fermi mixture 



]. We employ methods from conformal quantum 
field theory (CFT) to determine the asymptotic (long distance, low-energy) behaviour of 
correlation functions in the model from a finite size scaling analysis of the exact spectrum 
obtained by means of the Bethe ansatz. This approach gives the complete set of critical 
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exponents of the model as a function of the parameters in the microscopic Hamiltonian (see 



e.g. [6 



14. 



, ^ for applications to ID correlated electrons). As an application we compute the 
momentum distribution function of bosons and fermions in the atomic mixture as a function 
of their respective densities and the effective coupling constant. It should be emphasized 
that our results can be expected to describe the generic (universal) low-energy behaviour of 
atomic mixtures. Additional interactions - as long as they do not lead to a phase transition 
- will merely change the anomalous exponents but not the qualitative behaviour of the 
correlation functions. 

The ID Bose-Fermi mixture of N = Mf + M& particles with repulsive interaction (c > 0) 
on a line of length of length L subject to periodic boundary conditions is described by the 
integrable Hamiltonian |l2| 

N r) 2 

i=l 1 i<j 

Here Mf = Mf + Mj_ of the particles are fermions carrying spin a =f, \ and M& of them are 

e eigenstates of (JTJ are parametrized by the solutions of the Bethe 



bosons. The many-partic 
ansatz equations (BAE) 
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where e a (x) = (x + ia/2)/(x - ia/2) and M = N, M x = N - M T , M 2 = M b . The 
corresponding eigenvalue of (pQ) is i£ = X]^=i(<?j°' 1 ) 2 - I 11 ^ ne thermodynamic limit L — > oo with 
Mi/ L kept fixed, the root configurations {g^} of Eqs. (J2J) can be described by distribution 
functions pi which, as a consequence of (J2J) are solutions to |3] 

*(z)=pf +^£«®p i )(aO (3) 
i 

Here pf^ = (c/27r)<5j and are linear integral operators acting as {J. = dy) 

(K ij ®f)(x) = f k lJ (x-y)f(y) . (4) 
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The kernels of these integral operators are kij(x) = ai(x)5u-j\i — 02(^)^,1^,1 where 
27ra n (x) = 4n/(Ax 2 + n 2 ). The properties of the system are completely characterized by 
the densities = Mj/L = J.pW(x) of the components of the mixture (these relations de- 
termine the boundaries Qi of the above integral equations) and the dimensionless coupling 
strength 7 = Lc/N. For later use we also introduce the fraction of bosons in the system 
a = M b /N. 

Generically, i.e. for m; > 0, there are three modes of collective elementary excitations 
above the many-particle ground state of (0). Their dis per sion €i(k) is linear at low energies 
with different sound velocities Vi = dei/dk, i = 0,1,2 These quantities determine the 
finite size scaling behaviour of the ground state energy 



i ^ ' 

Physical excitations of the system are combinations of the elementary ones. Due to the 
interacting nature of the system the different modes are coupled and excitations in one of the 
modes shift the energies in the other ones. In general, this effect can be described in terms of 
generalized susceptibilities which may be determined in an experiment or numerically from 
studies of small systems 0. For the Bethe ansatz solvable models it is possible to describe 
the coupling of the modes in terms of the dressed charge matrix which in this case reads 

z ij = tij(Qi) ■ (6) 

The functions are given in terms of integral equations £ij(x) = 5ij + J2k(Kik ® £,kj){x)- 

Z determines the general form of the finite size corrections to the energies of low-lying 
excitations 

A£(AM, D) = QaM t (Z t ) ^ VZ^AM 

+ B T ZVZ T B + J2^k (iV+ + N-)) +o(^j . (7) 

Here, V = diag(t> , ^1, ^2) is a 3 x 3 matrix of the sound velocities, N k are non-negative 
integers, AM is a vector of integers denoting the change of Mj with respect to the ground 
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state for charged excitations. The Di are integers or half-odd integers according to 



D ~ (AM + AMi) /2 = AM T /2 mod 1 

Di ~ (AM + AM 2 ) /2 = AM//2 mod 1 (8) 
L> 2 ~ (AMi + AM 2 ) /2 = AM^/2 mod 1 

and enumerate finite momentum transfer processes: 
AP( AM, D) = ^ ( AM T . D + X) W ~ 

V k 

+ 2k FA Do + 2k Fd (D + D x ) + 2k B D j • (9) 

3 

Here k F ^ = irM a /L are the Fermi momenta of the fermion components, k B = irMb/L is the 



corresponding quantity for t 
In the framework of CFT 



re interacting bosons. 



3 the finite understood as that of 

a critical theory based on the product of three Virasoro algebras each having central charge 
C = 1 0, ^ | . Correlation functions of a general operator in the theory - characterized by 
the quantum numbers AMj and Di - will contain contributions from these three sectors. 
The simplest ones, analogues of primary fields in the CFT, have correlation functions (in 
Euclidean time t) 



exp (2iD kF j \X + 2i (D + Z?i) k F \X + 2% ( JT. D~ ) k B x 
(Mx, r)<M(0, 0)) = ^ — , — ~ t ' . , 2A - ' ~ ■ (10) 

The operators </>a are characterized by their scaling dimensions A^ in the chiral (left- and 
right moving) components of all three constituent theories. The latter are uniquely deter- 
mined from the finite size energies (J7J) and momenta © and form towers starting at 

2Af = ^Z kj D j ±±J2AM j (Z- 1 ) jl )j . (11) 

The asymptotic exponential decay of correlation functions in a large but finite system or 
at finite temperature T can be obtained from (|l()jl by conformal invariance. For exam- 
ple, at T > the denominators in (jl())l have to be replaced by by (v^t ± ix)~ 2A; 
(nT/v k sm tyT(t ± ix/v k )) 2Ak . 
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With ([11)1 the critical exponents which determine the long-distance asymptotics of any 
correlation function are known as soon as we have computed the dressed charge matrix 
(jSJ. To calculate the correlation functions of a given local operator O in the microscopic 
theory one needs to know its expansion in terms of the fields 0a of the CFT. Usually, 
this expansion is not known but O and 4>a have to generate the same set of selection rules 
in calculating the correlation function. This drastically reduces the number of possible 
terms in the expansion: As an example consider the bosonic Green's function Gb{x, r) = 
(^b{x, T)W b (Q, 0)): clearly generates a state with AMj = 1 which implies AMj = 1 in 
By (JHJ) the quantum numbers Dj are further restricted to integers: the uniform part 
of the Gb(x) ~ |cc | i/ 2 ^i. j s described by the operator with Dj = which allows to identify 
the TL parameter Kb [18] from (jlljl . The interactions lead to additional contributions to 
Gb oscillating with wave numbers k = 2kp a , 2kg, ■ ■ ■ 

For a comparison with experimental data one is often interested in Fourier transforms of 
the two-point correlation functions given above. The large distance behaviour of (jlO)) deter- 
mines the singularities of spectral functions near uj « ±Vk(k — k ) (see e.g. p). Quantities 
accessible in experiments with cold gases are the momentum distribution functions of the 
constituent particles. For the bosons this is the Fourier transform of the equal time Green's 
function Gb{x). From (fTU|) its singularities at wave numbers ho are then nb{k) ~ \k — k \ Ub 
near k « /c - The exponent z/& is the minimal value of 2 (Efc^fc + ) — 1 compatible 
with the quantum numbers AM and the selection rules for the D for the given ko, e.g. 

l/2K b = v h + 1 = | E fc (Ei(^"*)) a f o r k o = 0. 

Using the same procedure for the fermionic Green's function G a (x,T) = 
(^> a (x, r)\l/J.(0, 0)) we find that their asymptotic behaviour is determined by the confor- 
mal fields with AM = 1, AMi = AM 2 = 0, half-odd integers D , Di and integer D 2 for 
G T and AM = 1 = AM 1( AM 2 = 0, half-odd integers D 1 , D 2 and integer for G x . 
Again, the singularities of the fermions' distribution functions n a (k) follow from (fT0|) . Near 
ko — kpa = 0, ±2ks, ■ ■ ■ they are given by 

n c {k) ~ sign(A; — k )\k — k \ Uf , for k m k . (12) 

Vf is related to the dimensions (fTT|) for the quantum numbers AM and D as v b above. The 
Fermi distribution of non- interacting particles corresponds to Viikvn) = 0. 



In the following we consider two cases of particular relevance |ll|, [12J, 2 , namely (i) the 
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unpolarized case where Mf = Mj_ = Mf/2 and the ground state of the system is invariant 
under rotations in the spin index of the fermions and (ii) the fully polarized case where there 
is only one spin component of the fermions. 

The unpolarized gas. For Q\ = oo one obtains Mf = Mj. from (jHJ), i.e. with vanishing 
net magnetization. In this case the dressed charge matrix (jH} takes the form 



z = i 

2 



2Coo (£ 00 + ( i) 2( i ' 



(13) 





^ 2 Cio (do + Cn) 2 Cn J 

Here the Wiener-Hopf method has been used to determine Z\\ = l/y/2 and Coi = (oj(Qo 
Cij = Cij{Q2)- The functions Qj( x ) are given by 



Cij(x) = 8^ + R(x- y)( 0j (y) + R(x- y)dj(y) (14) 

JO 72 

with i?(x) = (l/7r) J °°dc(j e - '^/ 2 cos(cjx)/ cosh(u;/2). Using ()13)1 the scaling dimensions 
in (jllj) are independent on the remaining system parameters, i.e. the effective coupling 7 
and the bosonic fraction a. This a consequence of the £77(2) invariance of the system in this 
case. The mode ei(k) is the spinon mode of the unpolarized system, the CFT describing its 
low-energy properties is an [£77(2)] 1 Wess-Zumino-Witten model. 

Additional simplifications arise in the strong coupling limit 7 — > 00 (i.e. Qq — > 0) where 
Coo = 1) C10 = 0, C10 = Cn(0) — 1 = ol and Cn( a; ) is given by a scalar integral equation resulting 
from ()14j) . In Fig. Q we present results obtained from the numerical solution of these integral 
equations for the exponents which determine the singularities of the momentum distribution 
functions for bosons at k = and fermions at k = kp as a function of the bosonic fraction 
a for various values of 7. The exponents v^j at the other wave numbers are always larger 
than 1. Note that the system is in a different universality class for a = or 1. At a = 0, all 
particles are fermionic and the critical exponents are those of the ID Fermi gas jf|. Here the 
exponent Uf for the singularity at the Fermi point varies between and 1/8 as a function of 
7. On the other hand, the limit of i/& as a — > 1 gives exactly the exponent of the ID Bose 
gas with 5-interaction [lsj ]. 

The spin-polarized gas. Setting Q 2 = 00 in Q corresponds to Mj = 0. This case has 
been discussed recently in Ref. Q where some correlation functions have been computed 
numerically in the strong coupling limit. In this case, the finite size spectrum and the 




FIG. 1: Exponents characterizing the singularities in the bosonic (upper panel) and fermionic 
(lower panel) momentum distribution function for the unpolarized gas at k = and &f, respectively, 
as a function of the bosonic fraction a in the mixture for 7 = 0.2, 1.0, 5.0, 25.0, 00 (bottom to 
top). 

scaling dimensions are determined by two gapless modes. Again, the equations simplify in 
the strong coupling limit where all exponents can be given as a function of a directly, e.g. 
Vb{0) = a 2 /2 — a and Vf(kp) = o? — ct + 1/2 for the dominant singularities of the bosonic and 
fermionic momentum distribution functions, respectively (see Fig. E] for the 7 dependence). 
While the dependence of Uj, on a is similar to the one found in the unpolarized case, the strong 
coupling behaviour of i/f at small bosonic fraction is seen to be very different. Note that the 
singularity at kp + 2kg becomes very pronounced for sufficiently small a (yf = a 2 + a + 1/2 
at strong coupling). This feature of the fermionic distribution function is a direct signature 
of the interaction and should be observable in experiments. 

In summary we have used predictions from CFT on the finite size scaling of the low-energy 
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FIG. 2: Same as Fig. ^ for a mixture with polarized fermions. 



spectrum to compute the critical properties of a ID Bose-Fermi mixture for the integrable 
model In the generic case there are three linearly dispersing modes which determine 
the low-energy effective theory and the asymptotic behaviour of the correlation functions. 
Within this formalism we have related the critical exponents directly to the parameters 
describing the microscopic Hamiltonian, i.e. the coupling strength, the fraction of bosons 
and polarization of the fermions. The approach can be used to investigate the phase diagram 
of the ID mixture by identifying the order parameter with the slowest long-distance decay 



of its correlation functions (smallest exponent) 
is no instability leading to a phase transition 12 







111 ]. Within the integrable model there 



131 ] . The general expression (fTT|) . however, 



allows to obtain estimates on the exponents in a more general system based on numerical 
data on the spectrum of finite systems. 

This work has been supported by the Deutsche Forschungsgemeinschaft. 
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